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An  unlimited  source  of  messages  transmits  via  a  single  channel 
Poisson  wise  (  X  );  message  holding  times  are  random  (ExpOO). 
Conflict  and  mutual  -destruction  may  occur:  an  attempt  to  transmit 
when  channel  is  busy  destroys  both  incumbent  and  applicant;  both 
revert  to  a  limbo  state,  reapplying  independently  according  to 
Exp (  u  ).   Surprisingly,  this  process  may  reach  statistical 
equilibrium  with  a  simply  solution,  which  is  display 
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EXECUTIVE  SUMMARY 

Analysis  of  procedures  for  allocating  channels  to  randomly 
arriving  message  traffic  by  mathematical  methods  aids  in  making 
decisions.   This  paper  provides  an  example.   First,  we  provide 
some  background. 

There  are  situations  in  which  many  different  message  sources 
independently  compete  for  service  at  one,  or  a  finite  number  of 
channels  or  other  facilities.   If  one  message  is  on  the  channel, 
and  another  applies,  a  conflict  occurs,  and  both  are  destroyed, 
meaning  that  they  go  into  a  re-try  or  limbo  state,  from  which  they 
independently  attempt  to  access  the  channel  at  random  (exponentially 
distributed)  intervals;  of  course  retries  can  also  collide  and  be 
destroyed,  but  the  process  continues  indefinitely. 

The  usual  formulation  assumes  that  messages  appear  in  discrete 
packets;  then,  assuming  an  infinite  message  source,  and  that 
interrupted  packets  must  re-start,  the  number  in  limbo  will  eventually 
increase  indefinitely  --  the  process  is  unstable,  with  delays 
increasingly  to  infinity.   Only  by  allowing  for  takedown  or  defection 
in  the  model  can  stability  be  reached. 

This  report  generalizes  the  classical  model  slightly  to  consider 
a  long  stream  of  very  small  packets.   Interruption  is  still  possible, 
but  re-start  is  not  required.   It  is  shown  that  for  this  setup 
stability  may  be  achieved  provided  demand  rate  is  less  than  a  critical 
value.   The  probability  distribution  of  limbo  state  takes  on  a  simple 
form,  as  do  the  low  moments  (mean  and  variance) . 


A  RESOURCE  CONFLICT  RESOLUTION  PROBLEM 
FORMULATED  IN  CONTINUOUS  TIME 

D.P.  Gaver 
G.  Fayolle 


1 .   Introduction 

In  many  situations  involving  data  transmission  from  diverse 
sources   there  can  be  conflict  for  a  limited  number  of  channels 
or  other  facilities.   Uncoordinated  attempts  by  several  sources 
to  use  a  single  facility  can  result  in  "collision,"  the 
"destruction"  of  all  participants  in  the  collision,  meaning  the 
loss  of  the  transmission,  and  hence  the  need  for  re-transmission . 
An  important  problem  concerns  the  development  of  workable  proce- 
dures for  alleviating  the  conflict  and  corresponding  message 
delay  problems. 

Often  such  problems  are  viewed  as  occurring  in  discrete  time: 
slots  of  equal  length  occur  in  temporal  succession,  and  each  slot 
can  handle  just  one  packet  of  data  at  a  time,  if  two  or  more 
packets  try  to  use  the  same  slot  simultaneously,  a  collision 
occurs  that  somehow  must  be  resolved.   A  recent  paper  by  Fayolle, 
Flagolet,  and  Hofri  (1983),  hereafter  FFH ,  analyzes  a  stack 
protocol  for  handling  such  a  situation,  but  there  are  many  other 
proposals . 

This  report  is  concerned  with  some  simple  models  for  conten- 
tion for  a  single  facility  (channel),  and  for  contention  or 
conflict  resolution.   The  models  are  formulated  in  a  continuous- 
time  manner:   messages,  or  numbers  of  packets  constituting 
messages,  are  "long,"  meaning  that  they  occupy  many  consecutive 
slots  on  the  average  if  a  single  transmission  is  occurring. 


2 .   Model  1:   Poisson  Message  Source,  Single  Facility 

Permit  messages  to  arrive  at  a  single  facility  (e.g.,  bus  or 
satellite  link)  in  Poisson  manner  with  rate  A.   Service  times  are 
IIDExp(p).   When  a  message  arrives  it  either:   (i)  encounters  a 
free  facility  and  immediately  begins  transmission,  or  (ii)  inter- 
rupts ("collides  with,"  "destroys")  a  message  in  progress;  the 
result  is  that  both  messages  are  affected,  and  some  retransmis- 
sion becomes  necessary. 

In  what  follows  we  investigate  a  scheme  to  allow  the  messages 
to  retransmit  following  interruption  or  collision.   It  can  be 
called  a  stochastic  stacking  procedure.   In  a  general  sense  it 
is  patterned  after  the  algorithm  analyzed  by  FFH . 

2 . 1   Stochastic  Stacking  Model;  A  Single  Limbo  State 

Introduce  this  procedure:   whenever  a  collision  occurs,  each 
message  source  selects  a  delay  time  for  retransmission  indepen- 
dently from  the  distribution  Exp(v).   While  any  message  experi- 
ences the  latter  delay  it  will  be  said  to  occupy  a  limbo  state; 
the  number  of  messages  occupying  the  limbo  state  at  time  t  will 
be  denoted  by  X(t) .   Furthermore,  let  A(t)  denote  the  state  of 
the  facility  at  t:   A(t)  =  1  if  the  facility  is  occupied  with  a 
message  transmission,  while  A(t)  =  0  if  the  facility  is  idle  at 
time  t.   The  idea  of  the  delay  time  here  qualitatively  resembles 
the  randomization  scheme  studied  by  FFH. 

Apparently  tA(t)  ,  X(t)  ;  0  >_    t  j  is  a  Markov  chain  in  continu- 
ous time.   The  forward  Kolmogorov  equations  for  the  process 
can  be  written  in  terms  of: 


q.  (t) 


P(X(t)  =  j,A(t)  = 1}  , 


(2.1a) 


p.  (t)   =   P{X(t)=j,A(t)=0}  , 


(2.1b) 


The  probabilities  in  question  actually  depend  upon  initial  condi- 
tions, i.e.,  values  of  X(0)  and  A(0) ,  but  these  will  be  left 
implicit.   Consider  the  evolution  of  the  probabilities  as 
follows : 

p.(t  +  dt)  -  p  .  (t)  [l-Adt-jvdt]+q  .  „(t)Adt+q.  .  (t)v(j-Ddt 
3  3  D-2         ^3~1 


+  q  .  ( t)  udt  +  o(dt 


2.2a) 


q.(t+dt)  =  q . (t) [l-Adt-jvdt-udt]+p. (t) Adt 
3  3  3 


+  p.  , (t)v(j+l)dt  +  o(dt) 
j  +  1 


(2.2b) 


Now  subtract  p.(t)  (q.(t))  from  both  sides  of  (2.2a)   (2.2b), 
divide  by  dt  and  let  dt — ^ 0  to  obtain  the  formal  Kolmogorov 
equations : 


dt 


-  (  >.+vj)  p.  (t)  +Aq  ._2  (t)+v(  j-l)q  ,_  1(t)+uq  .  (t) 


(2.3a) 


dt 


-  (A  +  vj+y)q  .  (t)  +.-p  .  (t)+v  (  j  +  l)p.+1 


(t) 


(2.3b) 


A  few  words  of  explanation:   in  order  for  the  system  to  be  in 
state  (j,0)  at  time  t+dt  it  must  have  either:   (1)  been  in  that 


state  at  t,  i.e.,  a  moment  before,  and  have  experienced  no 
change,  or  (2)  been  in  state  (j-2,1)  a  moment  before,  and 
experienced  an  exogenous  (A-rate)  arrival;  the  latter  collides 
with  the  message  on  the  facility,  and  both  enter  the  limbo  state, 
so  the  state  changes  to  (j,0) ,  or  (3)  been  in  state  (j-1,1)  a 
moment  before  and  experienced  an  endogenous  (from  limbo  state) 
arrival,  or  (4)  been  in  state  (j,l)  and  experienced  departure 
of  the  message  on  the  channel.   This  explains  equation  (2.2a); 
otherwise,  in  order  that  the  system  be  in  state  (j-1,1)  at  time 
t+dt  it  must  have  either  (5)  been  in  that  state  at  t,  i.e.,  a 
moment  before  and  thus  experienced  no  change,  or  (6)  been  in 
state  (j,0)  and  experienced  an  exogenous  arrival  which  begins 
service  on  the  channel,  or  (7)  been  in  state  (j+1,0)  a  moment 
before  and  experienced  an  endogenous  arrival,  from  a  message 
that  leaves  the  limbo  state  and  entered  the  idle  facility.   This 
explains  expression  (2.2b). 

2 . 2   Long-Run  or  Steady  State  Conditions  and  Distributions 

To  look  for  the  conditions  allowing  a  stable  long-run 
distribution  of  {X(t),A(t)}  set  the  time-derivatives  to  zero 
in  (2.3a  and  b)  and  introduce  generating  functions  for  the  limit- 
ing probabilities;  by  definition 


P(z)   =   lim   7   p.(t)zD   =       p.zD,  (2.4a) 

t-»*  j  =  0   D  j  =  0   D 


Q(z)   =   lim  I      p.(t)zD   =    I      q.zD  ,  (2.4b) 

t-*oo  j=Q   J  j  =  0   J 


where  the  latter  limits  are  assumed  to  exist  at  least  when 
|z|  <_  1.   Performing  the  summations  leads  to  these  equations 
for  the  generating  functions: 

AP(z)  +  vzP'(z)   =   [Az2+u]Q(z)  +  vz2Q'(z)  (2.5a) 

AP(z)  +  vP'(z)   =   (A+y]Q(z)  +  vzQ'(z)  (2.5b) 

Now  multiply  (2.5b)  through  by  z,  subtract  from  (2.5a),  and  divide 
by  (1-z)  to  obtain 

AP(z)   =   (u-Az)Q(z)  .  (2.6) 

If  we  put  z  =  1,  this  expression  results  in 

A  [P(l)  +Q(1)  ]   =   yQ(l)   , 

but  since  we  assume  an  honest  limiting  probability  exists, 
i.e.,  that 


P(l)  +  Q(l)   =    I     [P{X=j,A=0)  +  PlX=j,A=.l}]   =   1  , 

j  =  0 


it  follows  that 


I   q.   =   Q(l)   =   -   =   p    (0  <  p  <  1)  .     (2.7a) 
j=0   3  y 


and  that 


I      p.   =   P(l)   =   1-p  (2.7b) 

j=0   3 


Further  information  results  by  differentiating  (2.6)  and  substi- 
tuting into  (2.5b);  the  result  is  the  simple  differential 
equation 


Q'  (z)      A  ( 1  +  z )  +  v      A  ,A (1  +  z)  +  v-| 
Q(z)         ,p_  _ 2z)       v   p  - 2Az 

A 


(2.8) 


J-  +   -(A  +v  +  £) t*- 

2v    v        2  u  -  2A  z 


which  can  easily  be  integrated  to  yield 


Zn^KTTT^       =      T^d-z)  +  ^U  +v  +^)£n(-^^-) 

Q(l)        2  v  2v         2     p  -  2  A  z 


or,  utilizing  the  value  of  Q(l)  , 


Q(„   =   P[e2v      (f^&>     2V    )  .  (2.9) 

1  -  2pz 


so  the  ergodicity  condition  immediately  appears  to  be  0  <  p  <  ^ 
It  then  follows  from  (2.6)  that 


P(z)   =   (i-z)Q(z)  ,  (2.10) 


and  hence  the  generating  function  of  total  system  occupancy 
(server  plus  limbo)  is 


H(z)   = 


E[zX+A] 


P (z)  +  zQ(z) 


'e2v        <r^>i      2v 


(2.11) 


0  <_  2p  <  1  . 

Differentiation  of    H(z)  at  z  =  1  generates  cumulant-like 
objects  that  can  be  converted  to  moments.   These  result  in 


E  [X  +  A] 


I^V1+2^  v">  ; 


(2.12) 


Var [X  + A] 


p(l  +p(3-p)^) 
(1  -2p)2 


0  <  2p  <  1 


(2.13) 


Apparently  both  the  mean  and  variance  of  total  system  occupancy 
decrease  monotonically  with  a  decrease  in  v   ,  the  mean  delay 
time  selected  by  (or  for)  any  colliding  message.   This  suggests 
that  the  best  control  policy  is  to  insist  that  interrupted 
or  destroyed  messages  should  immediately  try  again  for  trans- 
mission, if  one  is  guided  by  an  assessment  of  mean  system  delay 
time  by  Little's  Formula. 


2 . 3   Inversion  of  H(z)  When  Limbo  Delays  Are  Short 

Suppose  v  -»■  °°  in  (2.11),  signifying  retransmission  after 
a  negligible  delay  in  the  limbo  state.   Then  H(z)  approaches 

u  ,  .      ,  1-  2p  ,1/2  _  x  .     .  (2.14) 

H0(Z)   =   (1  -2p2}     ;   0  <  2p  <  1  . 

The  latter  resembles  the  ordinary  M/M/l  queue  generating  function, 
but  the  power  1/2  instead  of  1  is  noticeably  different.   Recog- 
nition that  the  generating  function  is  that  of  a  particular  nega- 
tive binomial  distribution  yields  the  explicit  formula 

r<:  +t)  1/2  j 

p{X+A=j}       =      f-(l-2p)  (2p)        ;       j    =    0,1,2,...  (2.15) 

ji    r<f) 

Suppose  for  one  moment  that  perfect  information  were  available, 
and  that  arriving  messages  could  be  queued  before  transmission 
on  the  system;  no  collisions  can  occur.   Then  the  probability 
distribution  of  the  total  number  in  the  system  is  well-known  to 
be  geometric, 

P{X+A=j)   =   (1-p)  pj  ,   0  <  p  <  1        (2.16) 


and 


E[X+A]   =   t-^—  (2.17) 

1-p 


It  then  appears  from  Little's  formula   that  the  ratio  of  long-run 
expected  total  delays  in  the  collision-prone  but  stacked,  and  the 


queued  systems  is  at  best 


'   P  0    <    p    <    1/2 


E [Delay    Stacked]  p  (?    ,R. 

E [Delay   Queued]  ]  l  ' 

1  +/O  +  20     +  oto    )       as      p    -    0    . 

which  clearly  reveals  the  advantage  obtainable  if  information 
can  somehow  reduce  or  remove  collision  frequency. 


3.   The  Busy  Signal  Problem. 

Consider  the  following  classical  problem,  historically  pre- 
ceding the  conflict  problem  previously  discussed,  but  of  interest 
in  its  own  right.   A  telephone  line  serves  a  number  of  customers. 
At  moments  of  a  Poisson  process  of  rate  A  customers  attempt  to 
initiate  calls  on  the  line;  if  the  line  is  free  it  is  captured  by 
a  caller  for  an  exponentially  distributed  time,  mean  y   .   If  a 
call  is  in  progress  when  another  call  arrives  the  newcomer  hears 
a  busy  signal,  hangs  up  and  tries  again  (retries)  after  an  exponen- 
tially distributed  time,  mean  v   ;  he  continues  to  re- try,  along 
with  others  who  have  experienced  busy  signals,  in  such  a  manner, 
i.e.  at  independent  exponential  (v)  intervals  until  he  accesses 
the  line  and  can  initiate,  and  eventually  complete,  his  call. 

The  above  familiar  setup  is  very  similar  to  the  conflict 
resolution  problem  just  addressed,  but  has  no  "collision"  or 
"destruction"  features.   Note,  however,  that  certain  proprietary 
telephone  systems  do  have  a  destruction  feature  for  low  priority 
calls.   Thus  the  U.S.  Dept.  of  Defense  AUTOVON  system  allows  high 
priority  calls  to  displace  ones  of  low  priority.   It  seems  reason- 
able to  model  this  latter  situation  using  a  limbo  or  retry  state 
much  as  we  did  the  previous  conflict  problem.   The  present  discus- 
sion models  only  the  single-priority  setup. 

3 . 1   Probabilities  for  A  Single  Limbo  or  Retry  State. 

Again  consider  the  vector  Markov  chain  {A(t)  ,X  (t)  ,  t  _>0  }  ,  where 
A(t)  =  1  if  the  line  is  occupied,  and  =  0  if  the  line  is  free, 
while  X(t)  is  the  number  in  the  limbo/retry  state.   Again 


q,(t)   =   P{X(t)  =  j,A(t)  =1}  (3.1, a) 


D 


Pj(t)   =   P{X(t)  =  j,A(t)  =0}  ;  (3.1,b) 


we  will  write  formal  Kolmogorov  equations  to  describe  the 
evolution  of  the  probabilities: 

Pj(t+dt)   =   Pj (t) [1-vjdt-Adt]  +  q. (t)ydt  +  o(dt)  (3. 2, a) 

q.  (t+dt)   =   q.  (t)  [1-vjdt-Adt-pdt]  +q._1(t)Adt 

+  p.  1(t) v(j+l)dt  +  p.(t)Adt  +  q.(t)vjdt  .     (3.2,b) 

The  usual  steps  yield  the  differential  equations 

dp. 

-^-     =   -(A  +vj)p_.  (t)  =  yq^(t)  (3. 3, a) 

dq 

-t£      =   -(A  +  y)q  (t)  +  Aq    (t)  +  v(j+l)p    (t)  +Ap  (t)  .    (3.3,b) 

J         j    ^-  J  J 

Notice  that  the  effect  of  one  possible  change,  i.e.  that  in  which 
a  retry  population  number  retries  and  again  gets  a  busy  signal, 
can  be  removed,  for  there  is  no  net  change  in  state. 

3 . 2   The  Long-Run  Distribution. 

Assume  now  that  a  long-run  distribution  occurs,  and  search 
for  the  necessary  conditions  and  the  distributional  form.   We 
must  solve  the  balance  equations  obtained  by  zeroing  the 
derivatives  in  (3.3) : 


(A+vj)p.  =   uq.  (3. 4, a) 


(A  +  p)q^   =   Aqj_]_  +  v(j+l)pj+1  +  Ap_.  .  (3.4,b) 


Introduce  generating  functions 


P(z)   =    I      z^p   ,   Q(z)   =    I      zDq   .  (3.5) 

j=0     J  j=0     J 


With  a  little  calculation  it  can  be  shown  that 

AP(z)  +  vzP'(z)   =   yQ(z)  (3. 6, a) 

(A+y)Q(z)   =   AzQ(z)  +  vP'(z)  +  AP(z)  (3.6,b) 

by  multiplying  (3. 4, a)  and  (3.4,b)  through  by  7?    and  summing, 
as  before;  the  primes  denote  z-dif ferentiation .   To  solve, 
multiply  (3.6,b)  by  z  and  subtract  from  (3. 6, a)  to  obtain 

A[P(z)  +  zQ(z) ]   =   uQ(z)  (3.7) 

after  division  by  1-z;  whence 

A[P(1)  +  Q(l)]   =       A   =   yQ(l) 


so 


lim  P{A(t)  =1}   =   Q(l)   =   A   =   p  ,  (3.8) 


the  probability  of  a  busy  line  in  the  long  run  (p  <  1) .   Next 
rewrite  (3.7)  as 

P(z)   =   (1/p  -  z)Q(z)  (3.9) 

and  differentiate, 

P'(z)   =   (1/p  -  z)Q'(z)  -  Q(z)  ;  (3.10) 

finally  substitute  for  P  and  P'  from  (3.9)  and  (3.10)  into  (3. 6, a! 
to  obtain  the  differential  equation 


Q*(z)   =   i/p  YVZ   Q(z)  (3.11) 


which  is  immediately  solved  to  yield 

,       1+A/v 
Q(z>   =   P(|  _  PnJ  (3.12) 

an  expression  for  P(z)  comes  from  (3.9),  and  finally 


,  .   „r  X(t)+A(t),      _,  .     ~  ,    >      ,1-  p N 1+A/v  , -,  ,_,> 
lim  E[z  ]   =   P(z)  +  zQ(z)   =   U ^— )    /  .  (3.13) 

t-^oo  L      pz 


It  follows  by  inspection  that  the  stationary  distribution  of 
system  occupancy,  including  the  channel  occupant  and  retry 
population,  is  now  negative  binomial: 


lim  P{X(t)+A(t)    -j>      =      jlnUA/t?'       Pjd-P>1  +  A/V    •  (3-14) 


t-*-°° 


Note  that  if  v  -*  °°  the  generating  function  of  (3.13)  tends  to 
( 1-p) / (1-pz) ,  that  of  the  geometric  distributionof  occupancy 
of  the  M/M/l  system  that  permits   queueing.   This  is  quite 
intuitive,  for  infinitely  frequent  retries  look  to  the  system — 
if  not  the  customer — exactly  as  if  arrivals  are  queued.   There 
is  a  decided  difference  between  (3.13),  or  even  (3.13)  with 
v  ->  °°,  and  the  corresponding  collision-destruction  model,  with 
generating  function  (2.11)  or  (2.1  ). 
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SOME  COMMENTS  ON  A  "RESOURCE  CONFLICT  RESOLUTION 
PROBLEM  FORMULATED  IN  CONTINUOUS  TIME" 

Alan  We  is  s 

AT&T  Bell  Laboratories 
Murray  HM,  New  Jersey  07974 

It  is  well  known  that  stoned  Aloha  schemes  are  strongly  unstable  when  there  are  an  infinite 

number  of  sources.  Yet  Gavex  and  Fayolle  [G-F],  in  the  paper  referred  to  in  the  title,  find  that,  for 

2p  <  1,  a  continuous  time  Aloha  scheme  is  stable.   In  this  note  we  try  to  explain  this  discrepancy 

via  some  elementary  calculations  which  illuminate  the  structure  of  the  model  analyzed  in  [G-F]. 

We  use  the  notation  of  [G-F]. 

Consider  X(t )  -  A(t),  the  total  number  of  messages  in  the  system  at  time  r .  We  shall  see  that 
this  process  is  approximately  a  birth-death  process  with  birth  rate  X  and  death  rate  p/2  when 
X  (r )  -  A  (t )  is  large.  Furthermore,  we  wQ  show  that  A  (t)  -1  about  half  the  time  when  X  -  A 
is  large,  and  we  will  develop  some  further  consequences  of  our  viewpoint. 

To  begin  with,  suppose  that  X(0)  -  K  is  large  and  that  A(0)  -0.  After  a  time  which  is 
distributed    exponentially    with    mean    —  =  — — ,   either   a    message    goes    from    hmbo    to 

transmission    (with    probability     — —  ~  1  —  — — )    or    a    new    message    arrives    (probability 

—  ~  ~rr- ).    Obviously,   the    most   likely   event    is   that   a    message    goes   from   limbo  to 

X  *  A  v         A  v 

transmission.      In     that     case,     after     another     exponemiaDy     distributed     time     with     mean 

— =   — — ,  either  another  message  from  limbo  collides  with  the  first  (probability 

X  —  (A  -  l)v  ^  ^l        A  v 

-— - =  1  -  •£- — ),     the     first     message     completes     transmission     (probability 

X  —  (A    —  l)v  —  ji  Av 

=   ■£-—)■  a  a   ncw   message  arrives  and   collides  with   the   first   (probability 


X  -  (A*  -  l)v  -  n        Av 


X  -  (A  -  Ijv  -  p.         Av 


=   -rr-).   In  any  case  we  see  that  in  a  (short)  time 


*  =  f 

A  v 


we  have  a  new  arrival  with  probability 


•2 


P  (kw)  *  YZ  *  *A/  , 
a  successful  transmission  with  probability 

p  (success)  *^=^, 
and  a  return  to  the  initial  condition  X  -K,A  ■  0  with  probability 

/>  (return)  =  1  -  XAr  -  ^  A/  . 

That  is,  the  process  is  approximately  birth-death  with  birth  rate  X  and  death  rate  y/2.  Furthermore. 

the  channel  is  idle  for  times  of  length  =  ■=—  and  is  busy  for  times  of  length  =  tt— ,  so  that  half 

Kv  Kv 

the  time  is  spent  busy  (this  is  another  way  of  saying  that  the  death  rate  is  ^).  This  clearly  shows 

2\ 
that  the  stability  condition  is  —  -  2p  <  1.  While  the  preceding  argument  is  not  rigorous,  it  is  not 

hard  to  tighten  up,  but  in  light  of  [G-F]  there  seems  to  be  no  point  in  doing  so. 

Let  us  carry  our  reasoning  a  bit  further.  We  have  (X(t)  +  A(t ))  is  approximately  a  constant- 
coefficient  birth-death  process  for  large  values  of  X  +  A ;  hence  in  steady  state,  we  should  have  as  K 
becomes  large, 

P{X(t)+A(t)-K)  =  C(2p)< 
for  some  C  >  0.  Using  the  notation  in  [G-F]  this  b 

A+*,-i*C(2p)*.  (1) 

Also,  since  the  process  spends  about  the  same  amount  of  time  in  state  (X  -  K ,  A  -  0)  as  in 
(X  -K  -1,  A  -1),  we  should  have 

Pk  ~  <li-\  •  (2) 

Together  these  two  equations  imply  that 

*/>r  -  VRk  ~  Mk  -1 . 

which  is  the  key  equation  (2j5)  in  [G-F].  Furthermore,  equation  (1)  is  the  same  as  their  equation 

2     ,z 

2 16  if  we  replace  C  by  -rrr  VI  -  p. 


Now  our  reasoning  was  valid  as  long  as  K  v  was  large.  Instead  of  assuming  that  K  was  large,  we 
could  as  well  have  supposed  that  v  was  large.  In  this  case  we  have 

1  2AT  -1 


*^ 


Kv       (K  -  l)v       K(F  -  l)v 


p  (successes)  ~ 


-  — Jt—  * 


<*  -  l)v       2  -  L 
K 


The  birth  rate,  of  course,  remains  X.  Then  by  a  basic  formula  of  birth-death  processes 

r    x, 
P(state  -AT)  -C  J[  -*■    forsomeC>0 

we  have 


/>(*(*) +A(r) -*)-/*  +^.t 

*    X 

-en  - 


2T 


Cp* 


2-1 


r  +  f 


'4 


vniri 


(2p)' 


which  is  the  equation  (2 IS)  in  [G-F]  for  the  choice  C  -  Vl  -  2p. 


We  can  also  analyze  the  effects  of  changing  the  model.  In  [G-F],  any  transmission,  no  matter 
how  short,  reduces  the  remaining  length  of  a  message.  In  fact,  we  saw  that  when  X(t )  is  large,  the 
system  takes  very  small  Writes"  very  quickly.  Now  a  realistic  system  may  have  some  overhead,  say 
a  time  6,  which  must  pass  before  a  new  transmission  may  start.  Or  there  may  be  some  basic  unit  to 
be  transmitted  such  as  a  packet  or  length  &,  and  if  the  transmission  is  interrupted  in  the  middle  of  a 
packet  then  the  entire  packet  must  be  retransmitted.  Or  there  may  be  a  time  delay  ft  through  the 
system  so  that  when  an  interruption  takes  place,  the  last  ft  of  the  message  must  be  retransmitted  to 
assure  its  safe  delivery.  We  model  any  of  these  cases  by  assuming  that  transmissions  which  are 
shorter  than  a  fixed  time  ft  are  ineffective. 


We  shall  now  demonstrate  that  this  system  is  strongly  unstable;  that  is 


P(hmX(t)  -»)  -1. 
•-  ■ 

Furthermore  we  shall  show  that 

P  (total  number  of  successful  transmissions  is  finite)  -  1; 
that  is,  we  will  show  that,  with  probability  1,  there  is  a  random  finite  time  T  such  that  no  successful 
transmissions  take  place  after  time  T.  To  this  end,  consider  the  probability  that  at  least  one 
message  which  attempts  to  transmit  during  a  time  when  X(i)  +  A(t)  -If  finds  no  uninterrupted 
interval  of  length  at  least  ft.  This  is  bounded  by  the  probability  that  in  a  time  interval  of 
exponentially  distributed  (X)  length,  a  Poisson  process  with  rate  K v  has  no  consecutive  events  more 
than  &  apart;  that  is, 

f  X      V        * 

P (no transmission)  >  2  (1  -*"*")y*1  1  ~  r    +x 


Kv  +  k 


(1  -*-*•) 


-l  -*:«-*"  +  o(*-*") 

That  is, 

P (transmission)  <  K\e-**>  -KK*-**') 
Hence  by  a  Bcrel-CanteHi  argument  we  may  conclude  that 

p  (infinite  number  of  transmissions)  »  0 
Since 

2  &*-"'  +  o(*-Mo  <  » 

for  all  positive  values  of  N  ,K  and  ft. 

Suppose  now  that  the  overhead  ft  b  a  random  variable  instead  of  a  constant.  Then  if 
P(t>  -  0)  *  0,  the  queue  is  strongly  unstable,  although  for  some  distributions  of  ft  an  infinite  number 
of  messages  may  get  through.  This  can  be  seen  as  follows:  work  is  done  on  a  message  only  after 


overhead  is  completely  serviced.  If  the  overhead  distribution  has  a  continuous  density  f(X)  near 
X  -  0,  then  successful  completion  of  overhead  occurs  approximately  as  a  Poisson  process  with  rate 
/(0)/2.  After  each  successful  overhead  completion,  a  message  gets  about  \I\X  (r)  time  to  transmit. 
Now  X{t)  will  dearly  be  on  the  order  of  Xf  for  large  t,  so  successful  message  transmissions  will 
occur  at  a  rate  of  about  /(0)/2X/ v  for  large  t .  The  analysis  of  other  types  of  overhead  distributions 
is  straightforward  and  is  not  defined  here. 

Let  us  now  suppose  that  message  lengths  are  not  exponentially  distributed,  but  have  instead  a 
distribution  function  F{X).   Then  it  is  not  hard  to  see  that  if  the  failure  rate  of  F(X)  (namely 


F'P) 


-)  remains  above  a  levd  p.  where  2X/ji<  1,  then  the  system  is  stable.    This  follows 


l-F(X) 

directly  from  our  original  analysis.  As  our  final  model,  consider  the  differences  between  Gaver  calls 
"resume,"  "res tart,"  and  "new"  disciplines.  "Resume'*  is  the  original  model,  where  each 
transmission,  no  matter  how  short,  reduces  the  remaining  length  of  a  message.  "Restart"  is  an  all- 
or-nothing  situation.  Either  the  entire  message  gets  transmitted ,  or  if  interrupted  it  must  be  entirely 
redone.  "New"  is  where  each  message  picks  a  fresh  length  at  the  start  of  each  transmission, 
independently  from  a  distribution  F(X).  That  is,  each  message  length  a  a  new  iixl.  random 
variable  in  each  transmission  attempt.  This  might  arise,  for  instance,  in  a  simulation  where  the 
programmer  does  not  keep  track  of  each  message  length  separately,  but  instead  picks  a  new  length 
each  time.  The  first  case  (resume)  has  already  been  discussed.  "Restart"  is  easily  seen  to  be 
strongly  unstable  by  the  same  sort  of  reasoning  we  used  for  analyzing  overhead,  and  the  reader  b 
invited  to  fill  in  the  details.  We  shall  now  show  that  when  F(X)  has  a  continuous  density  f(X) 
near  0,  "new"  b  stable  if  and  only  if  2X//(0)  <  1.  That  is,  the  three  cases  are  quite  different ,  and 
so  simulations  must  be  done  quite  carefully.  The  analysis  is  nearly  identical  to  the  "resume"  case: 
in  each  time  interval  of  length  A/  =  UK  v,  we  have  a  probability  of  abut  /(0)A//2  of  a  successful 
transmission.  In  fact  the  probability  of  a  successful  transmission  is  half  the  probability  that  a 
random  variable  with  distribution  function  F(X)  is  less  than  an  independent  exponentially 
distributed  variable  with  mean  UK  v,  so 


P(waxau*)  ~  j  f~  «-*•»  f(s)ds 
s=/(0)/2JTv. 
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